arXiv:astro-ph/0211409vl 18 Nov 2002 


Astronomy & Astrophysics manuscript no. melendez February 5, 2008 

(DOI: will be inserted by hand later) 


IRFM temperature calibrations for the Vilnius, 
Geneva, Rl(c) and DDO photometric systems * * 


Jorge Melendez^and Ivan Ramirez^ 


^ Seminario Permanente de Astronomfa y Ciencias Espaciales, Universidad Nacional Mayor de San Marcos. 
Ciudad Universitaria, Facultad de Ciencias Fisicas. Av. Venezuela s/n, Lima 1-Peru. 
e-mail: jorge@astro.iag.usp.br; ivan@f astmail. ca 
^ Institute de Matematica y Ciencias Alines, Jr. Ancash 536, Lima 1-Peru. 

Received 02 October 2002 / Accepted 07 November 2002 


Abstract. We have use d the in frared flux method (IRFM) temperatures of a large sample of late type dwarfs 
given by Alonso et al. (1996a) to calibrate empirically the relations Te// = /(colour, [Fe/H]) for the Vilnius, 
Geneva, RI(c) (Cousins) and DDO photometric systems. The resulting temperature scale and intrinsic colour- 
colour diagrams for these systems are also obtained. From this scale, the solar colours are derived and compared 
with t hose of the solar twin 18 Sco. Since our work is based on the same Te// and [Fe/H] values used by Alonso 
et al. (1996b) to calibrate other colours, we now have an homogeneous calibration for a large set of photometric 
systems. 
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1. Introduction 


Photometric calibrations of effective temperature can be 
regarded as fundamental tools to compare observations 
with theoretical models. They are needed to transform 
theoretical HR diagrams into colour-magnitude diagrams 
(e.g. VandenBerg 2000 ) as well as to perform the colour 
synth esis of stellar populations (e.g. Bruzual & Chariot 
|l993| ). Chemical abundance studies of stars, on the other 
hand, use the effective temperature (?///) as a key pa¬ 
rameter, so that reliable methods to calculate them are 
wi dely r equired (e.g Melen dez fc Barbuy ^002 ; Nissen et 
al. 2002; Ramirez & Cohen 2002). 

Since the first published photometric calibration of 


Teff (Popper 1959) many other papers devoted to the 
study of the colour relations between T///, [ Fe/H] and 
even logg ha ve bee n written (e.g. Code et al. 1976 , Bell 
& Gustaffson 1989 , Bessell et al. 1998 , Houdashelt et al. 
2000|). Most of them, unfortunately, considered only a rel¬ 


atively low number of stars, the problem being more crit¬ 
ical at the metal-poor end. To make this problem even 
worse we shall mention that these different calibrations 
used different sets of Te// and [Fe/H] values, thus yield¬ 
ing inhomogeneous empirical temperature scales. 

Alonso et al. (1996b; hereafter AAM96b) solved the 
problem partially by performing calibrations for the 
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* Based on data from the GCPD. 


BVRI(j), JHK(tcs) and uvby-/3 colours with a large sam¬ 
ple of cool dwarfs covering the ranges —3.0 <[Fe/H]< -1-0.5 
and 3500 K< Te// < 8000 K. The effective temperatures 
for the stars of their sample were computed by means of 
the infrared flux method (IRFM) whilst the correspond¬ 
ing me tallic ities were taken from the Cayrel de Strobel 
et al. ( 1992| ) catalog and partially from previously pub¬ 
lished photometric calibrations as described in Alonso et 
al. ( |l996a ; hereafter AAM96a). 


Even though it is possible to obtain the effective tem¬ 
perature of a star by applying the calibrations mentioned 
above, in some cases the only colour available is (B—V), 
which is very sensitive to metallicity, or colours that be¬ 
long to photometric systems that lack a homogeneous cali¬ 
bration. In this sense, the aim of this work is to extend the 
calibrations given in AAM96b to other colours for which 
a large amount of photometric data exist, like the Geneva 
system. 


In this paper, we present the calibrations for the 
Vilnius, Geneva, RI(c) and DDO photometric systems ob¬ 
tained with the sample in AAM96a and compare them to 
previously published calibrations. In Sect. 2 the sample 
and adopted photometry are described. The calibrations 
and the comparison with other works are presented in 
Sect. 3 leaving the discussion of the resulting temperature 
scale to Sect. 4. Conclusions are summarized in Sect. 5. 
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Colour F;(colour)/F;(B-V) 


(V-S) 

0.70 

(Y-Z) 

0.52 

(Z-V) 

0.29 

(B 2 -V 1 ) 

0.75 

(B 2 -G) 

1.01 

(V-R)(c) 

0.56 

(B-I)(C) 

0.69 

(V-I)(c) 

1.25 

C(42-45) 

0.23 

C(45-48) 

0.31 


Reference 


Straizys 1 
Straizys 1 
Straizys 1 
Cramer ( 
Cramer ( 

1995 

1995 

1995 

199t 

199£ 

) 

) 

) 

) 

Savage & Mathis 1 
Savage & Mathis ( 
Savage & Mathis ( 

1979) 

1979) 

1979) 

McClure 

McClure 

197£ 

197£ 

) 

) 


Table 1. Reddening ratios for late type dwarfs. 


2. The sample 


Although the atmospheric parameters T^ff and [Fe/H] 
of the sample used in this work are well described in 
AAM96a, we shall mention some extra details on these. 

All the stars of the sample are supposed to be late type 
dwarfs and subdwarfs. Nevertheless, a revision of their 


log 5 values in the Cayrel de Strobel et al. (|200l| ) catalog 
and their spectral types in SIMBAD database enabled us 
to identify some of these stars as giants and others as early 
type stars. We also found some variable and double stars. 

Several stars, apart from those belonging to the groups 
listed above, were excluded from our calibrations due to 
anomalies in their colours: G245-032, G009-016, G043- 
003, G014-039, G141-019, HD181007, HR8832, G190-015, 
G217-008 (anomalous Vilnius colours); HR66, G026-012 
(anomalous Geneva colours); vA 560 (anomalous Cousins 
colours); HD45282, HR4623, G154-021 (anomalous DDO 
colours). Stars with kinematical assignment of metallicity 
were also excluded as in AAM96b. 

The photometry adopted in this work was obtained 
from the General Catalogue of Photometric Data (GCPD, 
Mermilliod et al. 1997). For the sake of completeness we 
compiled all of the colours available for each star but, 
obviously, only a subset of these proved to be sensitive 
to the effective temperature, while being nearly indepen¬ 
dent of metallicity. Colours satisfying this condition are: 
(V—S) and (Y—V) in the Vilnius system; (B 2 —Vi) and 
(B 2 —G) in the Geneva system; (V—R)(c), (R—I)(c) and 
(V—I)(c) in the RI(c) system and C(42-45) and C(42-48) 
in the DDO system. The reddening corrections if(colour) 
were computed from the if(B—V) values given in AAM96a 
and the reddening ratios listed in Table |^. Unreddened 
(Y—V) = (Y—Z) -I- (Z—V) colours were obtained from in¬ 
trinsic (Y—Z) and (Z—V) colours. A similar procedure was 
employed to obtain unreddened C(42-48) colours. 


The residuals of the fits as a function of colour index 
and metallicity are often shown to check if systematic er¬ 
rors are introduced by the calibration formulae and also to 
illustrate their ranges of applicability in colour and metal¬ 
licity. 

3.1. Vilnius colours 

The Vilnius photometric system was originally designed to 
facilitate a bi-dimensional spectral classification for stars 
of any spectral type, luminosity class and reddening; the 
selection of medium-band filters being done in order to 
observe faint stars. 

For the colour index (V—S), which seems to be the 
best Teff indicator in the Vilnius system, we found, using 
data for 120 stars, the following fit: 

Oeff = 0.248 -b 1.336(V-S) - 0.354(V-S)2 

-0.029[Fe/H] - 0.008[Fe/H]^ , (1) 

for which the standard deviation a{deff) = 0.020. In terms 
oi Teff this amounts to 104 K. 

The application ranges of Eq. are: 

0.40 < (V-S) < 1.20 for - 0.5 < [Fe/H] < -bO.5 , 

0.40 < (V-S) < 0.70 for - 1.5 < [Fe/H] < -0.5 , 

0.45 < (V-S) < 0.70 for - 2.5 < [Fe/H] < -1.5 . 

The stars used to obtain Eq. (0) whose temperatures 
calculated from the calibration formula differ by more 
than 2cr from their IRFM temperatures are: G273-001, 
G021-006, G084-029 and HR6136. 

Figure shows the sample and the residuals of this fit 
as a function of (V-S) and [Fe/H]. There is a slight ten¬ 
dency towards lower temperatures for stars belonging to 
the [Fe/H] ^ —2 group but no other trends are apparent. 

The (Y—V) colours, temperatures and metallicities of 
134 stars allowed us to obtain the calibration formula 

9eff = 0.330 -b 0.981(Y-V) - 0.033[Fe/H] 

-0.009[Fe/H]^ , (2) 

with a{0eff) = 0.024 (126 K). This formula is applicable 
in the ranges: 

0.40 < (Y-V) < 1.05 for - 0.5 < [Fe/H] < -bO.5 , 

0.40 < (Y-V) < 0.65 for - 1.5 < [Fe/H] < -0.5 , 

0.45 < (Y-V) < 0.70 for - 2.5 < [Fe/H] < -1.5 . 


3. The calibrations 

Following the methodology of AAM96b, we have used the 
parameter Oeff = 5040/Teyy instead of Teff to perform 
the calibrations. In general, the multi-parametric relation 
6eff = /(colour, [Fe/H]) is a second order polynomial. 
Nevertheless, in order to improve our results, we have only 
considered those coefficients that are larger than 3 (t. 


The following stars depart more than 2a from the 
mean fit: G119-052, G217-008, G029-023, G021-006, 
G019-027, G014-024 and HR6752. 

Three stars with [Fe/H] < —2.5 were included in the 
sample used to calibrate this photometric system but since 
they showed a systematic tendency towards higher tem¬ 
peratures, we decided not to consider Eqs. ( 0 ) and § 
valid for the metal-poor end. 
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Fig. 1. Left: Teff vs (V—S) observed for the metallidty ranges —0.5 < [Fe/H] < +0.5 (filled circles), —1.5 < [Fe/H] < 
—0.5 (open circles), —2.5 < [Fe/H] < —1.5 (squares) and —3.5 < [Fe/H] < —2.5 (triangles). The curves corresponding 
to our calibration for [Fe/H] = 0 (solid line), [Fe/H] = — 1 (dotted line) and [Fe/H] = — 2 (dashed line) are also shown. 
Right: residuals of the fit (ATg// = T^ff — T™™) as a function of (V—S) (for the metallicity ranges indicated in the 
lower right section of the three upper panels) and [Fe/H] (bottom panel). 


A [Fe/H] 

ATe// for (V-S) 

ATe// for (Y-V) 

CO 

d 

1 

I 

0 

+ 

39 

45 

-1 —> -1.3 

16 

18 

-2 —> -2.3 

8 

9 


Table 2. Mean variation Are///A[Fe/H] in K per 0.3 
dex for Teff ~ 5000 K [(V-S) ~ (Y-V) cs 0.70 for solar 
metallicity] as a function of [Fe/H]. 


The high values of the residuals of these fits for the 
cooler stars (Teff ^ 4000 K) as shown by Figs. and 
^ could be associated mainly to rough estimates of the 
metallicity of this group. 

The gradients ATe///A(V—S) and ATe///A(Y—V) 
corresponding to our calibrations are both nearly inde¬ 
pendent of metallicity and vary from 70 K per 0.01 mag 
for Teff - 6000 K [(V-S) - 0.51, (Y-V) ~ 0.52] to 
25 K per 0.01 mag for Teff - 4000 K [(V-S) ~ 1.05, 
(Y-V) ^ 0.95]. As a consequence, an observational error 
of 0.01 mag in the photometry implies an error of about 
1.2% for the hottest stars and 0.6% for the cool end when 
using formulae and i)- 

On the other hand, the variations ATe/j/Ape/H] 
slightly depend on colour for both cases but the depen¬ 
dence of their values for a given Teff with [Fe/H] is simi¬ 
lar (see Table H). For bluer colours (greater temperatures), 
this quantity tends to increase, especially for solar metal- 
licities. 


For both (V-S) and (Y-V) the gradient 
ATe///A[Fe/H] approaches zero as [Fe/H] —> —2 

which shows that our formulae are reliable for halo stars 
even if their metallicities have been roughly estimated. In 
addition, this effect is independent of the colour values. 

3.2. Geneva colours 

A large quantity of stars have been observed with this 
photometric system in just two observatories making the 
homogeneity of the data excellent. We have found that 
approximately 60% of the stars with UBV colours have 
also Geneva photometry. Therefore, it will be very useful 
to have calibrations for this system. 

For the 0e//-(B2—Vi)-[Fe/H] relation we found the fol¬ 
lowing fit: 

9eff = 0.629 + 0.644(B2-Vi) + 0.065(B2-Vi)2 
-0.033(B2-Vi)[Fe/H] - 0.029[Fe/H] 

-0.010[Fe/H]^ , (3) 

obtained with data for 258 stars. The corresponding stan¬ 
dard deviation amounts to a{9eff) = 0.016 (90 K). 

From Fig. ^ it is clear that Eq. i) is applicable only 
in the following intervals: 

0.10 < (B 2 -V 1 ) < 1.05 for - 0.5 < [Fe/H] < +0.5 , 
0.10 < (B 2 -V 1 ) < 0.60 for - 1.5 < [Fe/H] < -0.5 , 
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Fig. 2. The same as in Fig. for Tg// versus (Y—V). 
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Fig. 3. 


The same as in 


Fig. 1^ for Teff versus (B 2 —Vi). 



400 p 


<1 


-400 




— ° 0 



-3 


-2 -1 

[Fe/H] 


0.20 < (B 2 -V 1 ) < 0.60 for 


-2.5 < [Fe/H] < -1.5 , 


0.20 < (B 2 -V 1 ) < 0.50 for 


- 3.5 < [Fe/H] < -2.5 . 


The stars that depart more than 2a from the mean 
fit are: G237-072, HR5901, HD3567, G251-054, G022-024, 
BD -1-71 31, HD111980, G126-062, G055-044, HD101177, 
G076-068, G227-037 and HR6844. 


It is also evident from this figure that the formula is par¬ 
ticularly reliable in the range 0.2 < (B 2 —Vi) < 0.6 where The residuals of this fit as shown by Fig. ^ suggest 
the star number density is high. a slight tendency of Eq. (^) to overestimate the effec- 
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tive temperatures of the hottest solar-metallicity stars and 
those of the [Fe/H] ~ —3 group. 

With the same stars used to obtain Eq. i) we found 
for the (B 2 —G) colour index the following formula: 


e^ff = 0.838+ 0.501(B2-G) 
-0.044[Fe/H] 


applicable in 

the 

ranges: 


-0.30 < 

(B 2 - 

-G) 

< 1.10 

for 

-0.25 < 

(B 2 - 

-G) 

< 0.40 

for 

-0.20 < 

(B 2 - 

-G) 

< 0.40 

for 

-0.20 < 

(B 2 - 

-G) 

< 0.30 

for 


0.030(B2-G)2 
0.010[Fe/H]^ , (4) 

- 0.5 < [Fe/H] < +0.5 , 

- 1.5 < [Fe/H] < -0.5 , 

- 2.5 < [Fe/H] < -1.5 , 

- 3.5 < ]Fe/H] < -2.5 . 


The standard deviation a{9eff) = 0.014 (86 K). 

The following stars depart more than 2cr from the 
mean fit: G021-022, HD3567, HD111980, HR7386, G076- 
068, G251-054, HD181007, HD101177, G227-037 and 
HR6844. 


The sample and residuals of this fit are both shown 
in Fig. Here we find that Eq. gives higher tem¬ 

peratures (by ~ 100 K) for solar-metallicity stars with 
Teff> 6400 K. 

There is only a slight dependence on [Fe/H] for the 
ATf.ff/A(B 2 —G) values, which vary from 50 K per 0.01 
mag at (B 2 —G) —0.2 to 15 K per 0.01 mag at 
(B 2 —G) ^ 0.8. In contrast with this behavior, the gra¬ 
dient ATe/z/Ape/H] is higher for the hottest stars (it 
could be as high as 130 K per 0.3 dex at (B 2 —G) ~ —0.2 
for solar metallicities) and lower for the coolest ones (it 
is only about 50 K per 0.3 dex for (B 2 —G) ~ 0.8 and 
[Fe/H] 0). Almost for any value of (B 2 —G), the gradi¬ 
ent ATe/y/Ape/H] ap proach es zero as [Fe/H] —> —2. 

Following Straizys ( 1995 , p. 372) we have also per¬ 
formed a calibration for the parameter t = (B 2 —G) — 
0.39(Bi—B 2 ), which is almost independent of [Fe/H] for 
F, G and early K stars. Data for 245 stars belonging to 
this last group allowed us to obtain the following fit: 


e^ff = 0.768 + 0.578t + 0.426t^ - 0.020[Fe/H] 

-0.007[Fe/H]^ , (5) 


with a{9eff) = 0.012 (75 K). 
to the ranges: 

-0.10 < t < 0.45 

for 

- 0.5 

-0.05 < t < 0.40 

for 

- 1.5 

+0.00 < t < 0.45 

for 

- 2.5 

+0.00 < t < 0.30 

for 

- 3.5 


This equation is restricted 

< [Fe/H] < +0.5 , 

< [Fe/H] < -0.5 , 

< [Fe/H] < -1.5 , 

< [Fe/H] < -2.5 . 


The following stars deviate more than 2a from 
the mean fit: G014-024, G021-022, G076-068, G048-039, 
HR159, G227-037, G126-062, HD3567, HD111980 and 
HD101177. 


Unlike formulae (||) and (^, the present calibration 
does not produce systematic effects on any metallicity or 
colour intervals as shown in Fig. ^ From this figure it 
is also clear that the range 0 < t < 0.4 provides good 
estimates of T/// almost for any [Fe/H] value. 

The calibration formula (H) introduces a shift of about 
50 K for an error of 0.01 mag in t if —0.1 < t < 0.2 for 
any value of the metallicity. This effect decreases mono- 
tonically as cooler stars are considered. On the other hand, 
the gradient ATeff/At vanishes for [Fe/H] ~ —1.5 inde¬ 
pendently of the t value and increases its absolute value up 
to 80 K per 0.01 mag for the hottest stars going towards 
both solar and very low metallicities. 

3.3. Cousins colours 

Since spectra of late type stars are dominated by emis¬ 
sion at long wavelengths, it is relatively easy to observe 
them with infrared filters. In addition, infrared colours are 
good temperature indicators because there are only slight 
blanketing effects in this region of the spectrum. 

For the (V—R)(c)-0e//-pe/H] relation we found: 

= 0.516 + 1.085(V-R)(c) - 0.203(V-R)^c) 

-0.120(V-R)(c)[Fe/H] + 0.025pe/H] , (6) 

using data for 97 stars. The standard deviation a{9eff) = 
0.016, which is equivalent to 91 K. 

Eq. (^ is applicable in the ranges 

0.20 < (V-R)(c) < 0.90 for - 0.5 < [Fe/H] < +0.5 , 

0.20 < (V-R)(c) < 0.50 for - 1.5 < [Fe/H] < -0.5 , 

0.25 < (V-R)(c) < 0.45 for - 2.5 < [Fe/H] < -1.5 ; 

which are the colour index and metallicity intervals cov¬ 
ered by the sample, as shown in Fig. Note that we have 
not considered the —3.5 < [Fe/H] < —2.5 interval because 
the only three stars with [Fe/H] < —2.5 actually have 
[Fe/H] ~ -2.6. 

The stars whose IRFM T^ff depart more than 2a from 
the mean fit are: G180-058, G021-022, HD157089, HR159 
and G251-054. 

Errors in the metallicity estimative are not so impor¬ 
tant when using Eq. ®) if we are dealing with E stars 
since ATe///Ape/H] < 20 K per 0.3 dex. The gradient 
itself is independent of [Fe/H] for any value of the colour 
index and reaches a maximum of 50 K per 0.3 dex for K 
dwarfs. 

On the other hand, the observational errors on 
(V—R)(c) do affect the temperature estimate for F stars 
since the variations ATe///A(V—R)(c) can be as high as 
90 K per 0.01 mag for (V—R)(c) ~ 0.22. For the coolest 
stars this effect can be reduced to 30 K per 0.01 mag. The 
metallicity dependence of this quantity is only slight. 

We have 104 stars with measured (R—I)(c) colours 
which satisfy the following formula: 

9eff = 0.422 + 1.557(R-I)(c) - 0.563(R-I)^c) 

-0.140(R-I)(c)[Fe/H] + 0.055[Fe/H] , (7) 
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Fig. 4. The same as in Fig. for Tg// versus (B 2 —G). 
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Fig. 5. The same as in Fig. for Teff versus t. 


and cover the application ranges: 

0.20 < (R-I)(c) < 0.80 for - 0.5 < [Fe/H] <+0.5 , 

0.30 < (R-I)(c) < 0.45 for - 1.5 < [Fe/H] < -0.5 , 

0.30 < (R-I)(c) < 0.50 for - 2.5 < [Fe/H] <-1.5 . 


The corresponding standard deviation a{6eff) = 0.013 
which amounts to 80 K. 


The stars whose temperatures obtained from this fit 
differ more by than 2cr with their IRFM temperatures are: 
G089-014, HD193901, HD74000, G026-012, HR159, G105- 
050, G088-027 and HD3567. 
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The relation between Tg// and (R—I)(c) is not affected 
by the star metallicity in the range 0.35 < (R—I)(c) < 0.45 
where the variation ATe///A[Fe/H] ~ 10 K per 0.3 dex, 
which means an error lower than 0.25% in the estimated 
Teff. For (R—I)(c) < 0.35 and (R—I)(c) > 0.45 this gra¬ 
dient increases to 30 K per 0.3 dex. 


Except for the coolest stars, the gradient 
ATg/z/A/R—I)(Q) depends strongly on [Fe/H], in¬ 
creasing its absolute value from 90 K per 0.01 mag for 
[Fe/H] 0 to 120 K per 0.01 mag for [Fe/H] ^ —2 
at (R—I)(c) ^ 0.3 for instance. For (R—I)(c) 0.4 an 

observational error of 0.01 mag would produce an error 
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of 70 K when using Eq. (Q) which means 1.3%, whilst for 
(R—I)(c) ~ 0.7 this reduces to 30 K (0.75%). 

The best temperature indicator for the Cousins sys¬ 
tem is the (V—I)(c) colour. Data for 97 stars satisfy the 
relation: 

= 0.483 -f 0.617(V-I)(c) - 0.072(V-I)^c) 

-0.066(V-I)(c)[Fe/H] -b 0.023[Fe/H] 







-0.008[Fe/H]^ , 


in the 

ranges: 






0.40 < 

(V- 

-I)(c) 

< 1.70 

for 

- 0.5 < 

[Fe/H] < 

-bO.5 , 

0.40 < 

(V- 

-I)(c) 

< 1.00 

for 

- 1.5 < 

[Fe/H] < 

-0.5 , 

0.55 < 

(V- 

-I)(c) 

< 0.90 

for 

-2.5 < 

[Fe/H] < 

-1.5 . 


The standard deviation corresponding to Eq. is 
aiOeff) = 0.011, which means the lowest dispersion ob¬ 
tained in this work (64 K). 

The following stars depart more than 2 ct from the 
mean fit: G021-022, HD3567, G105-050, HR159, G026- 
012 . 

The gradient ATe///A[Fe/H] depends strongly on 
colour and metallicity. Just to give an example, for 
(V—I)(c) ~ 0.7, where the star number density is high 
(Fig. ATe///A[Fe/H] vanishes for [Fe/H] ~ —1.5 but 
can be as high as 60 K per 0.3 dex for solar metallicities. 
If we go to (V—I)(c) ~ 0.4 the gradient now vanishes for 
solar metallicities and is higher than 100 K per 0.3 dex for 
[Fe/H] - -2. 

The ratio ATe///A(V—I)(c) monotonically decreases 
in absolute value from 60 K per 0.01 mag to 15 K per 






-1 

[Fe/H] 


0.01 mag depending on the mean value of (V—I)(c). For 
(V—I)(c) 0.7 this amounts to 35 K per 0.01 mag. 

We should mention that from Figs. 110 and ^ there 
is a slight tendency towards lower temperatures if they 
are obtained from the corresponding calibration formulae 
for the hottest solar-metallicity stars. Despite this fact, we 
have found that our calibrations for the Gousins colours 
do not introduce considerable systematic errors. 

3.4. DDO colours 

Although data in this system were available for only 89 
stars and even though DDO colours are very sensitive to 
the star metallicity, the need for an empirical DDO tem¬ 
perature scale encouraged us to perform the following cal¬ 
ibrations. 

For the G(42-45) colour we found the following fit: 

9eff = 0.492 -b 0.635[C(42-45)] - 0.067[G(42-45)]2 

-0.073[Fe/H] - 0.018[Fe/H]^ , (9) 

applicable in the ranges: 

0.40 < G(42-45) < 1.50 for - 0.5 < [Fe/H] < -bO.5 , 

0.40 < G(42-45) < 0.74 for - 1.5 < [Fe/H] < -0.5 , 

0.42 < G(42-45) < 0.76 for - 2.5 < [Fe/H] < -1.5 . 

Its standard deviation a{0eff) = 0.018 (103 K). 

Stars whose IRFM temperatures depart more than 
2cr from the fit are: HR5447, HR5568, HD111980 and 
HD140283. 
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C(42-45) 

Fig. 9. The same as in Fig. for Tg// versus C(42-45). 


Finally, for the C(42-48) colour, we obtained: 

= 0.189 + 0.413[C(42-48)] - 0.060[Fe/H] 

-0.015[Fe/H]^ , (10) 

whose application ranges are: 

1.30 < C(42-48) < 2.70 for - 0.5 < [Fe/H] < +0.5 , 

1.32 < C(42-48) < 1.82 for - 1.5 < [Fe/H] < -0.5 , 

1.32 < C(42-48) < 1.82 for - 2.5 < [Fe/H] < -1.5 . 

The standard deviation of Eq. amounts to 

a{6eff) = 0.015 (83 K) and only three stars depart more 
than 2cr from the fit: HR5568, HD111980 and HD140283. 

The general behavior of the gradients 
ATe///AC(42-45) and Arg///AC(42-48) is very 
similar. Their absolute values decrease for lower values 
of [Fe/H] and Tgff, and they both reach a minimum 
between 10 and 20 K per 0.01 mag for late K stars 
(C(42-45) > 1.2, C(42-48) > 2.2). For the hottest stars 
these quantities depend slightly on metallicity, and are 
about 40 K and 30 K per 0.01 mag for C(42-45) and 
C(42-48), respectively. 

The high values found for the mean variations 
ATe///A[Fe/H] corresponding to Eqs. d) and Jl^ con¬ 
firm the strong metallicity dependence of DDO colours, 
which is also evident from the values of the coefficients 
for [Fe/H] and ]Fe/H]^ in the calibration formulae. For 
Teff ^ 6200 K and solar metallicities, for example, 
ATeff/A[Fe/R] ~ 200 K and 160 K per 0.3 dex for 
C(42-45) and C(42-48), respectively. For [Fe/H] ~ —1 
these variations reduce to aproximately 70 K. For lower 
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temperatures ATey//A[Fe/H] decreases monotonically so 
that for Tg// ~ 5300 K and solar metallicity it is about 
140 K and 100 K per 0.3 dex for C(42-45) and C(42-48), 
respectively. 


3.5. Comparison with other calibrations 

In Figs. ^ to |l^ the present calibrations for [Fe/H] = 0 
and [Fe/H] = —1 (solid lines) are compared to previ¬ 
ously published theoretical and empirical calibrations as 
described below. 


3.5.1. Vilnius system 


The empirical work of Hauck & Kiinzli ( 1996| ; hereafter 
HK96) for the relation between Tg// and (V—S) (shown 
in Fig. O-a as a dashed line) is based on effective tem¬ 
peratures obtained by several authors using direct and 
semi-direct methods (including the IRFM). They fitted a 
linear function of (V—S) on 9f.ff with 66 stars, some of 
them out of the range covered in this work (namely, stars 
hotter than 7000 K). Since HK96 did not consider metal¬ 
licity effects on their calibration formula, we have assumed 
[Fe/H] = 0 for the comparison. 


For the hottest stars ((V—S) ~ 0.40) the difference in 
derived effective temperatures in our work and in HK96 
is higher than 150 K, for the coolest ones ((V—S) ~ 0.80) 
it is about —140 K (these differences amount to 2.2% and 
2.9 % respectively) whilst it vanishes for (V—S) ^ 0.50 so 
that the slopes are quite different. This can be attributed 
to the calibration formulae used in each work; the observed 
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C(42-48) 

Fig. 10. The same as in Fig. ^ for Tg// versus C(42-48). 


curvatures in our fits, for example, are due mainly to the 
(colour)^ term. 


In Fig. I^.a-b we also show the intrinsic colour indices 
determined for luminosity class V stars of various spectral 
types by averaging observed indices of unreddened and 
dereddened stars of the same MK spectral type (Straizys 


1995, p. 438). The calibration of Tg// in terms of spectral 
type is also provided in Straizys ( 1995| ). The agreement is 
good in the range 5800 K< Tg// < 7000 K for both (V—S) 
and (Y—V) while at cooler temperatures, the maximum 
difference in (V—S) is about 150 K (2.9%). For the cali¬ 
bration in (Y—V) ther e is a shift of aproximately —250 K 
for the Straizys ( 1995 ) work with respect to ours in the 
range 0.56 < (Y—V) < 0.92 but both slopes are rather 
similar. 


3.5.2. Geneva system 


The calibration given by Kobi & North (1990) is based 
on synthetic Geneva colours obtained from older versions 
of Kurucz models. It is compared to our calibration for 
the colour index (B 2 —Vi) in Fig. ^c. As usual, they cor¬ 
rected the synthetic colours using standard stars so that 
they could be able to reproduce the observed colours, es¬ 
pecially for cool stars. This correction was calculated only 
for [Fe/H] = 0 but they assumed that it was also valid 
for lower metallicities. Nevertheless, their published re¬ 
sults for [Fe/H] < 0 are ambiguous and have not been 
considered here. 


Some of the standard stars mentioned above were 
used by HK96 to calibrate this colour index and also 
the (B 2 —G), hence the better agreement between Kobi 
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& North (199C) and HK96 works for the coolest stars. In 


total, HK96 used more than 140 stars for the empirical 
calibrations shown in Fig. j^.c-d. The tendency of both 
calibrations with respect to ours are similar. They give 
higher temperatures than ours for the coolest stars whilst 
around Teff = 6800 K the difference vanishes. The maxi¬ 
mum difference amounts to 150 K for the (B 2 —Vi) colour 
index and 120 K for (B 2 —G). In both cases this maximum 
is achieved at Tg// ^ 4800 K and is equivalent to 3.0 and 
2.5 %, respectively. 

In a recent work, Blackwell & Lynas-Gray ( 1998[ ) used 
the IRFM temperatures of a large sample of ISO (Infrared 
Space Observatory) flux calibration stars and their Geneva 
colours (B 2 —Vi) corrected by interstellar extinction using 
Hipparcos parallaxes and an average interstellar extinc¬ 
tion law coupled with extinction maps to calibrate the 
relation Tg// = /(B 2 —Vi) where / is a second order poly¬ 
nomial. Their results are shown in Fig. ^c as crosses. The 
agreement between our work and theirs is excellent. We 
should also mention that the effect of m etallicity is con¬ 
sidered in Blackwell & Lynas-Gray ( 1998 ) but apparently 
valid only for [Fe/H] > —0.6. 


3.5.3. Rl(c) system 


Houdashelt et al. ( 2000| ) have obtained colour- 
temperature relations with improved MARGS stellar 


atmosphere models (Gustafsson et al. 1975; Bell et al. 
1976) after putting synthetic colours onto the obser¬ 


vational systems by means of a sample of field stars 
with IRFM temperatures determined and a model for 
Vega, whose colours fix the zero point corrections. These 
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(B,-G) 

Fig. 11. a-b. Comparison of our calibrations (solid lines) with those of Hauck & Kiinzli (|1996) (dashed line) and 


Straizys (1995, p. 438) (open circles) for [Fe/H] = 0; c-d. Comparison of Eqs. (pD an d (y) for [Fe/H] = 0 (solid lines) 
to other calibrations: Kobi & North ( 1990D (filled circles), Hauck & Kiinzli ( |1996| ) (dashed line) and Blackwell & 


Lynas-Gray (p^ ) (crosses). 
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Fig. 12. a-d. Comparison of our calibrations for [Fe/H] = 0, —1 (solid lines) with those of Houdashelt et al. ( 2000[ ) 
for [Fe/ H] = 0 (filled circles) and [Fe/H] = —1 (open circles), B essell et al. ( 199^ ) for [Fe/H] = 0 (triangles), Hauck & 
Kiinzli ( 1996[ ) for [Fe/H] = 0 (open squares) and Carney et al. ( 19941 ) for [Fe/H] = 0 (crosses) and [Fe/H] = — 1 (filled 
squares). 
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relations are shown in Figs. |T^.a-d as filled circles for 
[Fe/H] = 0 and open circles for [Fe/H] = —1. It is worth 
mentioning that the colour calibrations were derived 
for [Fe/H] = 0 stars so that their results for lower 
metallicities are less reliable. Their calibrations for solar 
metallicities seem to agree with ours for Tg// < 6000 K 
and for the whole range of temperatures in the case of 
(V—I)(c). For the hottest stars the (V—R)(c) colours 
derived by Houdashelt et al. ( 2000| ) are too red whilst the 
(R— I)(c) colours are too blue. The maximum difference 
with our work is about 150 K. 

In Fig. I^.b we show as open triangles the relation be¬ 
tween temperature and synthetic colours computed from 
spectra obtained from overshoot (for Tg// > 6000 K) 
and no-overshoot (for Te/f < 6000 K) models as given 
by Bessell et al. ( |1998 ). The selection of overshoot and 
no-overshoot models for the different ranges of tempera¬ 
tures alowed us to improve the agreement with our work 
although it showed to be very small. For effective temper¬ 
atures below 5000 K our empirical calibration is ver y close 
to this one but for the hottest stars Bessell et al. ( 1998 ) 
provide (V—I)(c) colours that are too red (the difference 
is about 0.03 mag). The agreement is also good, though it 
is not shown, for the (V—I)(c) colours of the late K dwarfs 
calculated on the basis of NMARCS models. 

We also show in Fig. |^.c the empirical calibration of 
HK96 (open squares), which, though being based on the 
effective temperatures of only 26 solar-metallicity stars, 
agrees perfectly well with the results of this work. 

In their survey of proper motion stars, Carney et al. 


(1994) derived the metallicity dependent relation between 


Teff and (R—I)(c) that is shown in Fig. ^c for [Fe/H] = 0 
(crosses) and Fig. ^d for [Fe/H] = — 1 (filled squares). 
Their temperature scale was based upon “a comparison of 
spectrophotometric scans of the Paschen continuum with 
surface flux distributions computed with the same model 
atmospheres used to calculate the synthetic spectra.” In 
general, their temperatures are lower than those calcu¬ 
lated from our formulae. 


3.5.4. DDO system 

Finally, in Fig. is shown the comparison between our 
calibrations for the D DO colours and those given by 
Tripicco & Bell 


and Claria et al. (|1994| ). The for¬ 
mer of these (which is shown in Fig. as filled and open 
circles for [Fe/H] = 0 and [Fe/H] = — 1, respectively) was 
based on MARCS models and used spectrophotometric 
scans for several stars to put the synthetic colours onto the 
observational system, thus allowing a model independent 
treatment. Tripicco & Bell ( |1991 ) argue that the effect of 
metallicity on DDO colours is so strong that “DDO sys¬ 
tem is primarily a measure of line, rather than continuum 
effects.” This is something that can be easily seen in Figs. 
P and and was already discussed in Sect. 3.4. For a given 


duces a decrease of 0.08 mag in C(4 2-45) and 0.09 mag 
in C (42-48) in the Tripicco & Bell ( |1991 ) work. For the 
present calibration, the decrease amounts to 0.1 mag in 
both cases. In this range there is a reasonably good agree¬ 
ment but large systematic differences start to appear as 
we go both to lower and higher temperature ranges. In 
general, the slopes are not similar but the effect of [Fe/H] 
on DDO colours is well reproduced. 

Spectroscopic observati ons o f held and open cluster 
stars allowed Claria et al. ( 1994 ) to calibrate the relation 
C(45-48) vs. C(42-45) for [Fe/H] ~ 0. This colour-colour 
diagram was then used to determine the mean values of 
DDO colours as a function of spectral type so that the 
mean DDO effective temperature-MK spectral type rela¬ 
tion provided them the Tg// :DDO colours calibration. The 
triangles in Fig. show that our temperatures are always 
lower than those derived by Claria et al (1994), the max¬ 
imum difference being about 200 K for K5 stars. There is 
a reasonable agreement around 5000 K. 


4. The empirical temperature scale 

Solving Eqs. (1) to (10) for the colour index as a function 
of Teff and [Fe/H] we obtained a set of intrinsic colours 
for late type dwarfs. They are listed in Tables || to ^ for 
various metallicities. Because of their usefulness, these re¬ 
sults have been used to plot the colour-colour diagrams 
shown in Figs. |IJ to |l^ Also shown in these hgures are 
the observed colours for stars with —0.5 < [Fe/H] < -1-0.5. 
The agreement is rather good since the observed disper¬ 
sion can be attributed primarily to the amplitude of the 
metallicity range considered. 

For the sun {Teff = 5777 K, [Fe/H] = 0), our temper¬ 
ature scale yields: (V—S)© = 0.547, (Y—V)© = 0.553, 
(Ba-Vi)© = 0.365, (Ba-G)© = 0.069, t© = 0.161, 
(V-R)®) = 0.352, (R-I)fc) = 0.328, (V-I)® ^ = 0.686, 
C(42-45 )q = 0.643 and C(42-48 )q = 1.655. These re¬ 
sults are in good agreement with the solar colours com¬ 
piled in HK96 and almost with every calibration men¬ 
tioned in Sect. 3.5 for it is clear from Figs. ^ to that 
the differences vanish for Tg// ~ 5800 K. It is also inter¬ 
esting to note that these solar colours are similar to the 
colours of the “solar twin” 18 Sco (Porto de Mello & da 
Silva |T9^: (V-S) = 0.54, (Y-V) = 0.57, (B 2 -V 1 ) = 
0.385, (B 2 -G) = 0.089, t = 0.17, (V-R)(c) = 0.353, 
(R-I)((.) = 0.335, (V-I)(c) = 0.686, C(42-45) = 0.651 
and C(42-48) = 1.674. 


Teff in the range 5200-6200 K, the effect of metallicity is 
such that a change from [Fe/H] = 0 to [Fe/H] = — 1 pro¬ 


5. Conclusions 

We have derived the relation between Teff, [Fe/H] and 
several colours for different photometric systems and have 
found good agreement with previously published calibra¬ 
tions. We have given detailed information about the ap¬ 
plication ranges and made several warnings on the use of 
our calibration formulae. The standard deviations amount 
from 64 K for (V—I)(c) to 126 K for (Y—V) and stars 
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C(42-45) C(42-48) 


Fig. 13. Comparison of our calibrations for [Fe/H] = 0, —1 (solid lines) with those given by Tripicco & Bell ( |1991 ) for 
[Fe/H] = 0 (filled circles) and [Fe/H] = —1 (open circles) and Claria et al. (1994) for [Fe/H] = 0 (triangles). 



(V-S) 

(Y-V) 

Teff 

[Fe/H] = 0 

-1 

-2 

[Fe/H] = 0 

-1 

-2 

3750 

1.205 

- 

- 

1.034 

- 

- 

4000 

1.049 

- 

- 

0.948 

- 

- 

4250 

0.932 

- 

- 

0.872 

- 

- 

4500 

0.839 

- 

- 

0.805 

- 

- 

4750 

0.763 

- 

- 

0.745 

- 

- 

5000 

0.698 

0.673 

0.667 

0.691 

0.667 

0.661 

5250 

0.642 

0.619 

0.613 

0.642 

0.618 

0.612 

5500 

0.594 

0.571 

0.566 

0.598 

0.573 

0.567 

5750 

0.551 

0.529 

0.524 

0.557 

0.533 

0.527 

6000 

0.513 

0.491 

0.486 

0.520 

0.495 

0.489 

6250 

0.479 

0.458 

0.453 

0.486 

0.461 

0.455 

6500 

0.448 

0.427 

- 

0.454 

0.430 

- 

6750 

0.420 

0.400 

- 

0.425 

- 

- 


Table 3. Intrinsic colours in the Vilnius system as a function of [Fe/H], 



(B 2 -V 1 ) 

(B 2 -G) 

t 

Teff 

[Fe/H] = 0 

-1 

-2 

[Fe/H] = 0 

-1 

-2 

[Fe/H] = 0 

-1 

-2 

3750 

1.008 

- 

- 

1.080 

- 

- 

- 

- 

- 

4000 

0.898 

- 

- 

0.890 

- 

- 

- 

- 

- 

4250 

0.800 

- 

- 

0.726 

- 

- 

- 

- 

- 

4500 

0.711 

- 

- 

0.583 

- 

- 

0.456 

- 

0.443 

4750 

0.631 

0.578 

0.555 

0.458 

0.386 

0.357 

0.393 

0.379 

0.380 

5000 

0.557 

0.507 

0.487 

0.347 

0.276 

0.247 

0.333 

0.318 

0.319 

5250 

0.490 

0.442 

0.424 

0.247 

0.178 

0.149 

0.276 

0.260 

0.261 

5500 

0.428 

0.382 

0.367 

0.158 

0.089 

0.061 

0.221 

0.204 

0.205 

5750 

0.370 

0.327 

0.314 

0.077 

0.009 

-0.019 

0.167 

0.149 

0.150 

6000 

0.317 

0.276 

0.265 

0.004 

-0.064 

-0.091 

0.115 

0.095 

0.097 

6250 

0.268 

0.229 

0.220 

-0.063 

-0.130 

-0.157 

0.063 

0.043 

0.044 

6500 

0.222 

0.185 

- 

-0.124 

-0.191 

- 

0.013 

-0.010 

- 

6750 

0.179 

0.144 

- 

-0.180 

-0.247 

- 

-0.038 

-0.062 

- 

7000 

0.139 

0.105 

- 

-0.232 

- 

- 

-0.089 

- 

- 


Table 4. The same as in Table ^ for the Geneva system. 


whose IRFM effective temperatures depart by more than 
2cr from the mean fit have always been listed. 

The solar colours derived in this work have proved to 
be similar to the colours of the solar twin 18 Sco so now 
these colours can be used to look for new solar twins. 

Along with AAM96b work, our calibrations provide an 
homogeneous (IRFM) temperature scale for cool dwarfs 


which can be used to explore the capability of models to 
reproduce the observations. A similar extension for giants 
is also necessary. 
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(V-R)(c) 

(R-i)(c) 

(V-I)(c) 


[Fe/H] = 0 

-1 

-2 

[Fe/H] = 0 

-1 

-2 

[Fe/H] = 0 

-1 

-2 

3750 

0.922 

- 

- 

- 

- 

- 

- 

- 

- 

4000 

0.808 

- 

- 

0.732 

- 

- 

1.534 

- 

- 

4250 

0.712 

- 

- 

0.638 

- 

- 

1.353 

- 

- 

4500 

0.631 

- 

- 

0.563 

- 

- 

1.201 

- 

- 

4750 

0.561 

- 

- 

0.501 

- 

- 

1.071 

- 

- 

5000 

0.500 

0.466 

0.439 

0.449 

0.443 

0.438 

0.958 

0.899 

0.879 

5250 

0.447 

0.419 

0.397 

0.405 

0.403 

0.402 

0.859 

0.814 

0.803 

5500 

0.399 

0.377 

0.360 

0.366 

0.369 

0.371 

0.772 

0.737 

0.735 

5750 

0.356 

0.339 

0.326 

0.332 

0.338 

0.343 

0.694 

0.669 

0.673 

6000 

0.317 

0.305 

0.296 

0.301 

0.311 

0.319 

0.624 

0.607 

0.617 

6250 

0.283 

0.274 

0.268 

0.274 

- 

- 

0.561 

0.551 

0.567 

6500 

0.251 

0.246 

- 

0.249 

- 

- 

0.503 

0.500 

- 

6750 

0.222 

0.220 

- 

0.227 

- 

- 

0.451 

0.453 

- 

7000 

0.195 

- 

- 

0.207 

- 

- 

0.403 

0.410 

- 


Table 5. The same as in Table || for the RI(c) system. 



C(42-45) 

C(42-48) 

Teff 

[Fe/H] = 0 

-1 

-2 

[Fe/H] = 0 

-1 

-2 

4000 

1.423 

- 

- 

2.593 

- 

- 

4250 

1.260 

- 

- 

2.414 

- 

- 

4500 

1.122 

- 

- 

2.254 

- 

- 

4750 

1.002 

- 

- 

2.112 

- 

- 

5000 

0.898 

0.792 

0.756 

1.983 

1.874 

1.838 

5250 

0.805 

0.702 

0.667 

1.867 

1.758 

1.722 

5500 

0.724 

0.623 

0.588 

1.761 

1.652 

1.616 

5750 

0.650 

0.551 

0.517 

1.665 

1.556 

1.519 

6000 

0.584 

0.486 

0.453 

1.576 

1.467 

1.431 

6250 

0.524 

0.428 

0.395 

1.495 

1.386 

1.350 

6500 

0.470 

- 

- 

1.420 

- 

- 

6750 

0.420 

- 

- 

1.350 

- 

- 


Table 6. The same as in Table ^ for the DDO system. 
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Fig. 14. (V—S) versus (B 2 —G), colour-colour diagram 
showing lines of equal Tg// and metallicity for [Fe/H] = 
0,-1 and —2. Crosses correspond to observed colours for 
stars with —0.5 < [Fe/H] < -1-0.5. 



Fig. 15. The same as in Fig. M for t versus C(42-45). 



C(42-48). 
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Fig. 17. The same as 
(B 2 -G). 
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